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Abstract 

In this paper we provide a step towards the understanding of the O(n) bulk 
operator algebra. By using a mixture of analytical and numerical methods, we 
compute (ratios of) structure constants, and analyse the logarithmic structure 
of the transfer matrix. We believe that the O(n) model for a generic value of 
n = q + q~ l ( i.e . for q not a root of unity) provides a toy model of a bulk 
logarithmic CFT that is considerably simpler than its counterparts at q a root 
of unity. 
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1 Introduction 

The most striking phenomenon in statistical physics is maybe universality in the 
vicinity of a critical point. As a statistical system undergoes a second order phase 
transition, the correlation length diverges, various quantities become related through 
scaling relations, and many power-law divergences occur, governed by universal 
critical exponents. While this phenomenon is well explained in all dimensions by 
the ideas underlying the renormalisation group approach, it is extremely difficult 
to compute the universal properties in general. At the critical point the system 
becomes scale invariant, which is usually promoted to conformal invariance, and 
the system exhibits fractal behaviour. The corresponding field theory is called a 
conformal field theory (CFT). 

Unlike in higher dimensions, the bidimensional conformal symmetry has infinitely 
many generators. This makes the CFT approach extremely powerful, and many 
critical exponents have been calculated exactly. This is the case for loop models 
such as the O(n) loop model and the Q-state Potts model, where the full spectrum 
has been known for some time [I], including the critical exponents for geometric 
objects such as the fractal dimension of interfaces. These models are particularly 
interesting as they contain many primary fields which are non-scalar (i.e. fields with 
a non zero conformal spin). Both the O(n) model and the Q-state Potts model can 
be written as (non-local) loop models, which are well known to exhibit logarithmic 
behaviour. For this reason, the underlying CFTs are poorly understood beyond 
their spectrum. 

A crucial part of the CFT data is given by its Operator Product Expansions 
(OPE), which are characterized by two main ingredients. The first one is the set of 
fusion rules, and the second one is the set of all structure constants, i.e. the three 
point functions between primary fields. Together these define the operator algebra, 
which has to be associative in a consistent field theory. For the O(n) model and Q- 
state Potts model the operator algebra is mostly unknown, which means that we do 
not have a complete understanding of how to compute correlation functions. Only a 
handful of structure constants have been obtained through an analytic continuation 
of the Liouville theory El E], but this method does not work for fields with a non¬ 
zero magnetic charge: a simple example where it fails is the probability for three 
given points to lie on the same loop. 
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Our objective is to describe the full operator algebra of these models, and in this 
paper we present a step in that direction. We focus on the Temperley-Lieb (TL) 
model, which is the loop model description of the Q-state Potts model. Moreover 
we restrict ourselves to generic values of the loop weight n (i.e. n = q + g -1 with q 
not a root a unity). 

The first tool we use is the Coulomb Gas description of the TL model. We extend 
the analysis of crossing symmetry of four-point functions to the case of non-scalar 
primary fields. While a subset of the structure constants can be extracted this way, 
charge neutrality proves to be a strong limitation on the set of structure constants 
for non-scalar fields that can be accessed with this approach. Therefore we turn 
to a more powerful approach inspired by the bootstrap of Liouvillc theory [4j. It 
is a well known fact of CFT that the decoupling of a null vector yields a differ¬ 
ential equation for correlation functions. In conjunction with crossing symmetry, 
it is possible to compute ratios of structure constants. Using an extension of the 
bootstrap method to the case of non-scalar fields, we make predictions for many 
ratios of structure constants. An interesting remark is the appearance of logarithms 
in some four-point functions, as expected from such non-diagonalisable CFTs. We 
also compute numerically these structure constants on the lattice, and we find an 
excellent agreement with our analytical calculation. 

The logarithmic nature of the CFT raises a fundamental question regarding 
the bootstrap approach : do null-vectors decouple in the scaling theory of the O(n) 
model? Indeed, while in a unitary CFT the null vectors - i. e. states with a vanishing 
norm - always decouple, logarithmic CFTs numn are known to have a non- 
dehnite inner product, allowing states of zero norm to have a non-vanishing inner 
product with other states. This is typically what happens in non-trivial Jordan cells 
of a logarithmic CFT. These questions can be also be addressed on the lattice: for 
the TL model the Jordan cells can already be observed in finite size on the lattice 
model (see 0 CEDI CEB El H3]), and this provides a powerful tool to analyse the 
decoupling of null-vectors. 

The paper is organised as follows. In Sec. [2j we recall the definition of lattice 
loop models, and review the operator content of the corresponding compact boson 
CFT (with a background charge) in the scaling limit [131 H], In Sec. [3j we generalise 
the Coulomb-Gas PUDS] and conformal bootstrap mug approaches to determine 
a family of OPE coefficients in this CFT. In Sec. [4] we come back to the lattice to 
study the exact Jordan blocks appearing in the spectrum of the transfer matrix (or 
Hamiltonian) of the loop models, and then discuss the associated structures in the 
CFT. Moreover, we give an example of analytic calculation of the indecomposability 
parameter j3 in one of these Jordan blocks. Finally, we present some open questions 
and perspectives in Sec. [5] In the Appendix, some technical details of our analytical 
calculations are given, and our numerical method to compute OPE coefficients is 
outlined. 
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2 The O (n) and Temperley-Lieb loop models 

2.1 Lattice models 

The 0 (n) loop model on the honeycomb lattice [14] consists in configurations of 
non-intersecting closed polygons. A configuration C is given the Boltzmann weight 

W Q{n) (C) = K^ n N , (2.1) 

where i(C) is the total length of polygons in C, N c (C ) is the number of contractible 
loops in C, and N nc (C ) is the number of non-contractible loops in C (see Fig. [Tj). 

For fixed n,n both in the interval [—2,2], the O(n) loop model has two distinct 
critical regimes: (i) when K is at the critical value K c = 1 /\/2 + \/2 — n, for 
example critical polymers (n — 0), critical Ising domain walls (n = 1), XY model 
(n = 2); (ii) when K > K C: the model is in the universality class of the Q-state 
Potts model, with n = \[Q. In this paper, we will mainly focus on the universality 
class ii, since it has a simpler realisation, namely the Temperley-Lieb model. 

The Temperley-Lieb model is defined on the square lattice, and consists in com¬ 
pletely packed configurations of non-intersecting closed polygons, as shown in Fig. [2] 
A configuration C is given the Boltzmann weight 

Wtl(C) = n Nc{c) , (2.2) 

with the same notations as above. For completeness, we recall the relation [18] 
between the Potts and TL models, in the case of a planar domain, and for n = n. 
Consider the Q-state Potts model on the lattice jC, with spin variables {s^} and 
Boltzmann weights exp(J Yl(ij) The partition function can be graphically 

expanded using the identity exp(J5 SiiS ..) = 1 + vS Si}S ., where v = e J — 1. This gives 
the Fortuin-Kasteleyn partition function 

ry _ \ ^ y^)#connected components ..^occupied edges /o 

A Potts / J W " i 

clusters 

where the sum is over all possible subgraphs (or “cluster configurations”) of the 
lattice C. Each cluster configuration can in turn be associated uniquely to a loop 
configuration of the TL model on the medial lattice Af, i.e. the lattice connecting 
the midpoints of adjacent edges of £: see Fig. [2j Using the Euler relation, one gets 

= Q N - m X V (v/VQ) K - (C, VQ N ‘ (C> , (2.4) 

loop config. C 

where N s (jC) is the number of sites of £, N e (C) is the number of edges of C not 
crossed by the loops, and N^{C) is the number of closed loops. 

2.2 Scaling limit: compact boson CFT 

The emergence of a compact boson in the scaling limit is best described for the 
honeycomb O(n) model on a torus. Take a configuration C of the O (n) model, and 
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Figure 1: An example configuration of the 0(n) model on the honeycomb lattice 
embedded on a torus. 



(b) 



Figure 2: (a) An example configuration of the Temperley-Lieb model on the torus, 
(b) The same TL configuration, with the corresponding Fortnin-Kasteleyn clusters 
and dual clusters (the Potts spins live on the full dots). 
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give an orientation to each loop, independently of the others. Then, define discrete 
height variables {4>j} at the centers of faces of the honeycomb lattice, so that the 
oriented loops are their contour lines, with a step ±7r across each contour line: see 

Fig. [3j 

Note that, for a contractible loop oriented anti-clockwise (resp. clockwise), the 
algebraic number of turns to the left is +6 (resp. —6). Hence, the factor n Nc ^ 


in (2.2) can be distributed locally by assigning each left (resp. right) turn a phase 
factor e* x (resp. e _ * x ), with n = 2cos6y. 

In the scaling limit, it can be argued that the coarse-grained analog of <pj is 
a free scalar field 4>(r). Moreover, the free scalar field 0(r) has to be compact. 
Indeed, a given cycle of the torus may be crossed by a non-zero flux 2m of arrows 
(with m G Z/2 for honeycomb O(n) and m € Z for TL), which results in a defect 
5(j) = 2nm for the height variable. The effective model for both critical regimes of 
the O (n) model and for the TL model is thus a compactified boson 


w - s 


d 2 r (V0) 2 



for honeycomb O(n), 
for TL. 


(2.5) 


The renormalised value of the coupling constant g can be fixed in terms of n using 
e.g. exact Bethe Ansatz determination of one critical exponent [18]. It is found to 
be 


n = —2 cos tt g , 


1 < g < 2 for K = K c , 

0<g<l for K > K c or in the TL model. 


( 2 . 6 ) 


For non-contractible loops, the local phase factors defined above cancel each 
other, so one must introduce an additional topological factor related to the height 
defects 5(p and 8'<j> along the two cycles of the torus. Simple geometric arguments 
show that the partition function on the torus is given by [lj 

Z 0(n) = J2 [ S^Km [ D $ e ~ m COS t 7re 0 ( 2m ) A ’ ( 2 - 7 ) 

m,m / GZ/2 

S(f)=27vm' 

where the symbol A stands for the greatest common divider, and 

n — 2 cos7reo . (2.8) 

For the TL model, we consider a torus consisting of L x M lattice steps, with L 
and M both even. Then the defects must be of the form 5(p = 2nm, with m G Z. 
Using the same line of argument as for the honeycomb O(n) model, one gets the 
partition function 


Z TL- I 

m.m'EZ 


5(f)=2'Km 
5' 4>=2Trm' 


[.D0] e cos(7re 0 m A m!). 


(2.9) 


In both cases, the insertion of the topological factor changes the dependence of 
Z on the system size, and hence it affects the central charge, which becomes 


fie 2 

c=l-^ 


( 2 . 10 ) 
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Figure 3: An oriented loop configuration and associated height values on the dual 
lattice. Heights are given in units of n. 


To make contact with the Coulomb-gas notations pa uni, we introduce the back¬ 
ground charge —2 a 0 = eo/s/g, and the central charge is given by c = 1 — 24 ckq. 

2.3 Vertex operators 

To discuss the operator content of the CFT associated to the loop model, we shall 
use the geometry of the plane, with complex coordinates (z,z). The height field 0 
can be expanded on holomorphic and anti-holomorphic modes, and one defines the 
dual field 6 as: 


</>(-, z) = <p(z) + <p(z ), 0(z, z) = <p(z) - <p(z). (2.11) 

A chiral vertex operator V a , with conformal dimension h(a), is given by 

V a (z) = : exp [iy/4g a ip{z)\ : h(a) = a 2 — 2acto , (2-12) 

where : • • • : denotes normal ordering. A general vertex operator is of the form 

V a ,a(z, z) = : e*V^[av>«+a?(*)] • (2.13) 

Let us recall a heuristic argument relating vertex operators to “spin-wave” (electric) 
and “vortex” (magnetic) excitations of the height variable 0. First, if a — a — a e i, 
the vertex operator has the form V e i = : exp(i v /4^ ct e i0) : and is single-valued only 
if ^/Aga e \ is a multiple of the inverse compactification radius : 


(y. e \ 


k 

Ws 1 


k e 2Z for honeycomb O(n), 
fcgZ for TL. 


(2.14) 


A second interesting case is when a = 2cko — a = a mag . Then the vertex operator 
has the form 


V mag = : exp {iy/Ag a 0 0) x exp[iy/Ig (a mag - a 0 )9] : 


(2.15) 
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where the first factor cancels the background charge, and the second factor creates 
a defect of height Sep = 87r(o; ma g — a 0 )/\/4 g, which must be a (half-)integer multiple 
of 2n\ 


Q^mag «0 


m-y/g j rn G Z/2 for honeycomb O(n), 
me Z for TL. 


More generally, we introduce the notation for vertex operator charges 


mJg e 
ttme = tto---h 


2 '2 Jg' 

and the corresponding conformal dimensions 


e 

ttme = «0 H-X-h 


2^’ 






4^ 


hme hrne T me . 


(2.16) 


(2.17) 


(2.18) 


The limit n —> n corresponds to e 0 —>■ 1 — <?, which gives 2a 0 —* y fg — 1/y/g. Thus 
if one sets ck_ = yj~g and cx + = — 1/y fg, one gets 


2a 0 ——)■ ct_|_ T ol _ , ck_|_ck_ — —1, 

n—>n 

1 1 

«me » v(l - m)a_ + -(1 - e)a+ , 
n— yn Z Z 


(2.19) 

( 2 . 20 ) 


which is, for m, e integers, the Coulomb-gas parameterisation of the Kac table. 


2.4 Operator spectrum 

Let us go back to the toroidal geometry, and consider the transfer matrix (or Hamil¬ 
tonian) propagating along one cycle of the torus. The transfer matrix conserves the 
total flux of arrows in the time direction, which we denote 2m. We distinguish two 
types of sectors: 


1. The sector rn — 0 corresponds to the situation where no loop cycles in the 
transfer direction, ft contains the ground state, and the allowed excitations in 
this sector are of the electric type, 


OL — (X — OL o e , e — 6q T k , 


k G 2Z for honeycomb O (n), 
fcgZ for TL. 


( 2 . 21 ) 


and we denote the corresponding operator as Vo e . Note that, in the limit 
n —> n, the charge ao,e 0 +fc a i,i+fc> and we obtain the scalar degenerate fields 
Voe —» Vi,i +fc = 'hip+fc with k G 2Z (resp. k G Z) for the honeycomb O(n) 
model (resp. for the TL model). 


2. A sector with rn ^ 0 corresponds to 2m strings propagating in the time di¬ 
rection. The lowest energy state in this sector is the “purely magnetic” state 
with a = 2a 0 — a = a m0 , and the allowed excitations are of the form 


(cp cr) (cK me , ol 


m G —, e G Z/m 
m G Z + , e G Z/m 


for honeycomb 0(n) 
for TL 


( 2 . 22 ) 


—m,ej ? 



and we denote the corresponding operator as W me - This operator has con¬ 
formal spin h me — h_ mte = —me. Note that, for m > 1, “fractional” electric 
charges, e = p/m with p integer, are allowed. 

To summarise, for a transfer matrix with periodic boundary conditions, the states 
present in the scaling theory are the Vo e ’s in the zero-string sector, the W me ’s in the 
2m-string sector, and their descendants under the action of the Virasoro generators 
{T n , L n j. 

3 Conformal bootstrap for mixed electric/magnetic 
operators 

Throughout Sec. [3] we will restrict to the case n = n. For the sake of completeness, 
and to fix notations, we first review the Coulomb-gas approach [ T5l |16 | for the case 
of scalar electric operators, and present a generalisation to mixed electric/magnetic 
operators, also with the help of functional relations in the fashion of [4]. 

3.1 Coulomb-gas approach for purely electric operators 

3.1.1 Four-point function 

We denote V a (z,z ) =: e l vUa 0 (>>z) : a generic electric operator. Recall that V Q and 
V- 2 ao -a have the same conformal dimensions h = h = a 2 — 2a^a, and they represent 
the same operator <&h,h hi the CFT. We consider the four-point function 

C(z,z) = (V ai ( 0 )V a 2 (z, z)V m (l)V a 4 (oo)) , (3.1) 

and we assume that the charges {cp} satisfy the neutrality condition after insertion 
of (p — 1 ) screening charges V Q+ : 

ay + a 2 + 03 + 04 + (p — l)a + = 2ckq . (3-2) 

Under these conditions, the four-point function can be written 

C(z, z) = ( zz) 2aia 2 [{z - 1 )(z- 1)] 2 “ 2 “3 x G ( Z) z) , (3.3) 

where G(z, z) is a bilinear combination of p independent conformal blocks 

p p _ 

G(z,z) = EE Xk,k Fk(z |ay, a2, 03, CI4) T (^|cvi, «2, CI3, 0:4). ( 3 - 4 ) 

k= 1 k= 1 

3.1.2 Conformal blocks 

The conformal block ^.(^lay, a 2 , 03 , CI 4 ) is dehned as the contour integral over the 
positions {vj} of the (p — 1 ) charges 

^■(-l{a;}) = / II dv i II ~ 1 ) 6 ( U J ~ Z T x II ^ _ V ^ P ’ ( 3 ‘ 5 ) 

JCk j =1 j =1 i<j<i<p-i 
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V a2 (z) V a3 (l) V a2 (z) V a3 (l) 


(a) 


(b) 


Figure 4: Fusion diagram corresponding to: (a) the conformal block T k \ (b) the 
conformal block Tp. 


where 

a = 2a + ai , b = 2a + a3 , c = 2a + a.2 , d = 2a + a4 , (3-6) 

and p = a+, and the integration contour C k is shown in Fig. 11 1 . The behaviour of 
J-fc at £ —» 0 is 

z 2aia2 T k (z) ~ N k x Z -M«i)-M« 2 )+Kft) (3.7) 

z—>• 0 


where 6 k = «i + «2 + (k — l)a+, and the normalisation factor N k is given in (A.l). 
Thus, the conformal block T k corresponds to the fusion diagram shown in Fig. |4 ^l, 
and the term (k, k) in (3.4) corresponds to the fusions 


V Q1 x V Q2 —* V Pkih <- V Q3 x V< 


« 4 


(3.8) 


3.1.3 Dual basis 

Another basis of conformal blocks J r p{z\ai 1 a 2 , a 4 ) is given by the same integrals 

as in (3.5), but with the integration contour C k replaced by Cg (see Fig. 11 3 ), and 
they correspond to the fusion diagram of Fig. [ijy with the charge fit = ot 2 + ot^ + 
(t — l)a + in the fusion channel. The correlation function may be written 


G(z,z) = £jf (/ X ( (z|{ ai }) X(*'l{«.}) 


(3.9) 


LI 


The change of basis is given by the matrix A and its inverse A: 


Fk(z) = 22 A keFt{z ), kFt(z) = 22 MkkFk{z) , (3.10) 

i =1 k =1 

with the matrix elements of A and A given in ( A.2| -A.3). The expansion coefficients 
of G in (3.4) and ( |3.9D are related by 

Xt.i = 22 Am A kl X k ,k ) X k .i = 22 A lk , (3-11) 

k,k ijt 


where we have used the fact that the coefficients A k g and A kk are real. 


Hk 
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3.1.4 Monodromy conditions 

Under analytic continuation of C(z,z) as z goes around zero, the term (k,k) in the 
decomposition (3.3 3.4) acquires a phase e 2l7rAtfe U where [see (3.7) 

Vkk = KPk) ~ HPk) • ( 3 - 12 ) 

Note that p kk is simply the conformal spin of the vertex operator Vg k p, in the fusion 
channel. Similarly, under analytic continuation of C(z,z) as z goes around one, the 
phase for the term (£,£) in (3.9) is e 2l7rMtf , where 

fiti= h (Pi)~ h (Pi)> ( 3 - 13 ) 

which is, again, the conformal spin of the operator in the dual fusion channel. 

We require G(z, z) to be monodromy invariant: this selects the terms in the 
expansions ( 3.4[ ) and (3.9), such that p kk and are integers. There are in general 
several possibilities, depending on the values of the external charges {a*}, but we 
consider here the generic solution, where the selection rules are: 


k = k 


and 


= £. 


(3.14) 


In other words, we must have X kk = 0 for all k ^ k, and X^j = 0 for all £ ^ i. This 
yields the “diagonal form” 


G(z,z) = j2x t 


(3.15) 


k =1 


and the homogeneous linear system 


vc/?, ^44^ = 1 


(3.16) 


k =1 


The system (3.16) consists in p(p — 1) linearly related equations, for p unknowns 
{Ad,..., X p }. It can be shown that (3.16) is equivalent to the subsystem obtained 
by taking £ E {1,... ,p — 1} and £ = p : 


V £ G {1,..., p .1} , A k £ A kp X k — 0, 


(3.17) 


k =1 


which has the simple solution X k oc A pk /A kp . To get a meaningful set of coefficients, 
we need to work with the normalised conformal blocks: 


G(z,z) = ^S k 


k =1 


k Mz) 


S k =\N t \ 2 xX k 


(3.18) 


After some simple manipulations described in Appendix [A| we obtain 

ft=nMw+iw 7(1+ ‘ +i ' ,)7(1+c+ " ) 

3=0 

p—k—1 

x n ^iu+^p] 


7 [2 + a + c + (k - 2 + j)p\ 

7(1 + b + jph(l + d + jp) 


3 =0 


'y[2 + b + d+(p-k-l+ j)p\ ’ 


(3.19) 
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where we have used the notation 'j(x) = r(x)/T(l — x). Similarly, in the basis {J^}, 
we have the expansion 


G(z,z) = j>t 


t= 1 


Np 


(3.20) 


where Sp is obtained by changing k <->• i and a b in (3.19). 


3.1.5 Restricted case 

To lighten the notation, we define V rs ■— V Qrs . Let us consider a particular case 


of (3.1-3.2): 


C(z, z) = (V a (0)V 2ao - a (z, z)V lp ( l)V lp (oo)), (3.21) 

where a is a generic charge. Note that, after proper normalisation V a = V 2 Q0 _ Q , but 


we use the above convention to ensure the neutrality condition (3.2). The term k 

(3.22) 


in (3.18) corresponds to the fusion channel: 

V Q x V2 q 0 -o —> Vp k Vi p x Vi p , 


where /3 k = 2a 0 — 07 , 2 * 1-1 ■ The coefficients S k then read 

c c(p)i i\ TT re 1 i\ 1 7(1 + a + jp)7(l + a> + 3P) 

S k = S k (a,a)=|| 7 [(j + l)p] - 7[(fc + j)p] - 

./lY.n, 1 7 2 (1-(P-J-1)P) 

^ II T[(l UPj rr) 111 ' 1 \ 1 > 

7[2-(p + fc-j-l)p] 

where a = 2a + a, and a' = 2a + (2ao — a). 

3.1.6 Extraction of OPE coefficients 

The expansion coefficients S k are related to OPE constants by 
S k CX C(V a , Va, Vi,2fc-l) X C(Vl,2 fe -l, V lp , Vip) . 


(3.23) 


Using the result (3.23), and after a few simple steps described in Appendix |A| 
obtain the compact form 


(3.24) 
we 



(3.25) 


where a = 2a + a , and a' = 2a + (2a 0 — a). Similarly, the dual coefficients for (3.21) 
give the OPE constants: 


Sp OC C (Vip, V Q , Vq._|_(p+ 1—2£)a+/2) • 


(3.26) 
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In this case, the OPE constants are determined up to an overall constant (indepen¬ 
dent of q): 


C 2 {V lp: V a ,V a+qa+ ) cx 
e-2 

UrtU + Dp] 7(1 + a' + jp) ^[1 + a + (p - £ - 1 - j)p\ 
j =0 

p—£—l 

x II 7 2 [(i + !)p] l^ + a + jp) lil + a'+ {£-2-j)p], 
j =0 


(3.27) 


where q G { — 2^, 
with V 12 ' 


., ^~ 2 ~} j and f = 2±I — g. An important example is the fusion 


V 12 X V a —> V a ± a+ /2 ■ 


(3.28) 


Using (3.27), we get 


C ,2 (Vl2,Va,Vq-a + /2) _ 7(1 + a')7(l + « ~ P) 

U 2 (Vl 2 , Vq, V q+q+ / 2 ) 7(1 + 0 ) 7(1 + a ' — p) 


3.1.7 Comparison with the DOZZ formula 

As it was noted in M, some of the structure constants of CFTs with c < 1 
associated to statistical models satisfy the same functional relations as those of the 
Liouvillc theory defined by the Lagrangian 

£ = i ( W) 2 + ^e 2 ^, (3.30) 

with the coupling constant b continued to an imaginary value: b = i/3. In this 
context, the central charge reads: 

c = 1 — 6(/3 — 1//3) 2 - (3.31) 


The solution of the functional relations (conformal bootstrap) for the Liouville the¬ 
ory is given by [T^ |2U[ fT?t |gj 


C c (V ai ,V a2 ,V a3 ) = (3.32) 

congt x TW ~ ~ <4i)1pW - /3- 1 - ELi «») 

[ni’.iT„(/?-2 aj )T42/3-/?-i-2a j )] 1/2 

where aA = ag + ay — a*,, and the multiplicative constant is simply determined 
by Cc(V Q ,V Q , 1) = 1. To define the function Yp(x), it is convenient to introduce 
Q = j 3 + /3 _1 . On the interval 0 < x < Q, the function Y p(x) is given by 


T/j(*) 



(Q /2 - x) 2 e 4 


sinh 2 (<5/2 — x)| 
sinh y sinli ^ 


(3.33) 
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and the function is fully determined for x G R by the relation 


T M^) 


(3.34) 


This solution is usually called the Dorn-Otto-Zamolodchikov-Zamolodchikov 


(DOZZ) formula. It can be checked explicitly that the OPE coefficients (3.25 


3.29), 


found from the integral representation of the conformal blocks, coincide with (3.32). 


3.2 Coulomb-gas approach in the presence of magnetic op¬ 

erators 

We shall generalise slightly the Coulomb-gas approach described above to the case 
of vertex operators V aj(5 with a a. 

3.2.1 Four-point function 

We consider the correlation function 

(Poi,ai (0 )Vq 2 ,02 ^)l^a3,a:3 (1)1^04,04 ( l -^- 1 )) j 

subject to the neutrality conditions 

Or + Q?2 T a 3 T OL4 + (p — l)a_|_ = 2 q?o j 

ol\ + 0-2 + «3 T 04 + (p — 1 )n+ = 2ao ■ 

We can then write 

C(Z, Z) = Z? aia *z mi " 2 ( Z - 1 ) 4 “ 2 “3^ _ 1 )4a 2 a 3 x Qfo ) (3 37 ) 


(3.35) 


(3.36) 


with 


p p _ 

G(*,z) = EE X k j. Tk{z |q;i, «2, «3, a 4 ) J^(^|ai, a 2 , «3, a 4 ), (3.38) 

fe=l k =1 


where the TVs are the same conformal blocks as in (3.4), but the charges in the 
holomorphic and anti-holomorphic sectors can be different. In the dual basis we 
write 


p p ~ ~ __ 

G(z,z) = EE Xiji F t (z\a u ot 2 i 03 , 04 ) J~^z\cx\, 01 2 , 03 , 04 ). 

*=i f=i 

We introduce the notations 


a = 2a + ai, b = 2a + a 3 , c = 2a + a^ , d = 20+04, 
a = 2a + ci 1 , b = 2 a + a 3 , c = 2a + a:2 , d = 2 o + a 4 . 


(3.39) 


(3.40) 
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3.2.2 Monodromy conditions 


The term k, k in (3.38) corresponds to the fusions 


Vai,ai X V a2! a 2 ^ ^ ^03 ,«3 ^ "^04,04 


with the charges in the fusion channel: 


/3k — cci + 0 L 2 + (k — 1 )o+ , 
13^ — ot\ T ot 2 T {k — l)a+ . 


(3.41) 


(3.42) 


Under analytic continuation of C(z, z ) as z goes around zero, this term gets a phase 
factor e 2m/ife >U where 


Mk,k = [h((3k) ~ h(ai ) - h(a 2 )\ - [h(/%) ~ H a 1 ) ~ h(a 2 )\ 
= spin(V /3fc ^J - spm(V aiiai ) - spin(V a2 , a2 ). 


(3.43) 


Similarly, under analytic continuation of C(z,z) as 2 goes around one, the term £,£ 
in (3.39) gets a factor e 2l7r ^, where 

Me,l = ) - spin(V Q2 ,a 2 ) - spin(V a3 ,a 3 ) • (3.44) 

In general, to construct a monodromy invariant four-point function, one must select 
the terms such that /i fc j, and Jlgg are integers. 


3.2.3 Restricted case 

From now on, let us focus on the particular case: 

C(z, z) = {V a ^(0)V 2ao - a ,2ao-a{z, z) Vi p (l)Vi p (oo)), (3.45) 

where we require that 

(i) the operator V a ,a has integer conformal spin; 

(ii) the difference (a — a) G Zct_/2. 

Note that these two conditions are satisfied by the electric operators {Vyi+fc} and 
the 2m-string excitations {W me }- We then have 

a + c = a + c = 2 (p — 1), b = b = d = d = b p =( 1— p)p , (3.46) 

and the condition (ii) translates into 

a — a = c — cGZ. (3.47) 

In this context, we can again construct a monodromy invariant by selecting only the 
terms k — k and £ — £. We get a “diagonal” decomposition 

p 

G{z,z ) = Y, x k Mz\M)Mz\{ai}), (3-48) 

k =1 
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with coefficients X k subject to the linear system 


V l± A ki{{ui}) X k = 0 . 


(3.49) 


k =1 


The above condition (ii) ensures that this system still admits a solution^ which has 
the form ~ 


Xk x A^a ,})' 

Introducing the coefficients S k on normalised conformal blocks: 


G(z,Z) = £> 


T k (z\{ai}) T k (z\{ai}) 


x 


k= i 


iV fc ({«,}) N k {{ai}) ’ 


we obtain from (3.19): 


S k = \jS k p \a, a') Sl P \a,a ’), 


(3.50) 


(3.51) 


(3.52) 


where 


a = 2a + a , a' = 2a + (2a 0 — a ), a = 2a + a , a' = 2a + (2a 0 — ct), (3.53) 

and S^ > \a,a') is defined in (3.23). 


3.2.4 Extraction of OPE coefficients 

The above expansion coefficients give access to the OPE constants through the 
relation 

s k (X C(V a , & , V a , 5 , Vi, 2fc _i) X C(V li2k -i, V lp , V lp ). (3.54) 

Using the results for purely electric operators, we get 


(3.55) 



where C(V a , V a , Vi^fc-i) is the purely electric coefficient given in (3.25). Similarly, 
using the dual expansion we obtain 


(3.56) 



where q G { — ..., and C (Vi p , V a , V a +qa + ) was given in (3.27). In particu¬ 

lar, for p = 2, we get 


C 2 (V 12 .V«, a .V a . ¥ia ^) 


/y( 1 + a') 7(1 + a')7(l + a - p)y( 1 + a - p) 

c 2 (v 12 ,v a , a ,v„ += ± a+ ^) \j 

7(1 + 0)7(1 + 0)7(1 + a' - ,0)7(1 + a' - p) ’ 


(3.57) 


If (ii) is relaxed, then (3.50) is a solution of (3.49) only for l < t. 
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Like for purely electric operators, these coefficients are simply related to the DOZZ 
formula ( 3.32 ). C(V q . 1i q, 1 , V a2 ,a 2;^«3,<S3) Gc 0 Fa 21 F23)dx (Fq^ j Fa 2 j Fa3 )• How¬ 

ever, note that this simple form of the OPE coefficients (3.55-3.57) as geometric 
means of the holomorphic and anti-holomorphic parts is only valid under the as¬ 
sumptions (i) and (ii) of Sec. 3.2.3 For instance, the OPE coefficient of purely mag¬ 


netic operators (7(W 10 , W 10 , W 10 ) is expected to be finite, whereas the formula (3.32) 
with charges a* = a 10 gives an infinite coefficient. 


3.3 Functional relation on OPE coefficients 

3.3.1 General setting 

Throughout this section we move away from the O(n) loop model and consider a 
generic OFT with a central charge of the form 

(1 — g ) 2 2 

c = 1 - 6 v = 1 - 2Aal (3.58) 

with g irrational. We the following assumptions about this CFT: 

(i) There are no degeneracies in the spectrum. 

(ii) The spectrum contains a scalar field $12 (z, z ) with conformal dimension h\ 2 = 
h \2 = (3 — 2 g)/Ag. This field is degenerate at level 2, namely the descendents 

X12 — (L -2 — gL\) $12 and Xi ,2 = {L -2 — d^i) $12 (3.59) 

have zero norm. 

(iii) Finally the last assumption we make is that both X12 and X12 decouple from 
the theory. Note that, in unitary CFTs, this last assumption follows from the 
previous ones, but not in logarithmic CFTs, where a null vector can still have 
non-zero overlap with its logarithmic partner. 

It is well known that the decoupling of X12 yields a differential equation for any 
correlation function that involves the field $ 12 . In particular it yields a functional 
equation for the structure constants of scalar fields [3]. In this section we extend 
this result to the case of primary fields } L with arbitrary spins. The method is 
essentially the same as in ^j, and relies on locality and crossing symmetry of the 
four-point functions. 


3.3.2 Four-point function involving $12 

Let us consider a four-point function of primary fields 


G(z,z) = (A? hl) - hl (0)^ 12 (z,z)^ h3 - h3 (l)^ h4 - h4 (oo )). (3.60) 

The decoupling of x.vi implies the following (holomorphic) differential equation: 


gd 2 z 


1 1 

“ 3 - 7 

z z — 1 


d z + ( — - 


z~l 


-\- h 2 — ht 




G{z,z)= 0. (3.61) 
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Likewise, the decoupling of X 12 yields a similar antiholomorphic PDE, with hi —>■ hi, 
z — y z, d z —y d s . In order to write down the solutions of these differential equations, 
it is convenient to parametrise the conformal dimensions a la Coulomb gas, namely 

h = a (a — 2a 0 ), h = a(a — 2a 0 ) (3.62) 

We stress that this is purely a convenient parametrisation of the conformal dimen¬ 
sions (h,h), and that we are in no way working with a Coulomb gas CFT. Note 
that the two values a and a' = 2 «o — ol lead to the same conformal dimension 
h(a) = —aoi. As a consequence, there are four ways to parametrise a field with 
given left and right conformal dimensions h = h(a), h = h{a) : 

(a, a), (a 7 , a), (a, a 7 ), (a 7 ,o; 7 ). (3.63) 

In the following we will write $ aj6; as a short-hand for &h(a),h(a)- At this point we 
put no constraint on the spectrum, but it will turn out that the decoupling of X 12 
and X 12 severely constrains the admissible values of (a, a). 


3.3.3 Conformal blocks 


Let us focus for now on the holomorphic conformal blocks, which are given by the 
two independent solutions of the PDE (3.61). The Coulomb gas parametrisation of 
the conformal dimensions is very convenient to write down explicit expressions for 


the conformal blocks. For this purpose we adopt the same notations as in Sec. 3.2 


a — 2a + oi\ , b — 2 crcr 3 , d — 2 cr_|_cr 4 , 

a 7 = 2 a + aq , b' = 2 a + a 7 3 , d' = 2a + a' A . 


(3.64) 


We will restrict ourselves to the generic situation in which the exponents at the 
singularities of the PDE do not differ by an integer. This means a — a 7 , b — b' and 
d — d' 2Z. This ensures that the conformal blocks have no logarithm in their 
expansion)^ 

Several basis of conformal blocks can be considered. I 11 the following we will use 
two such bases. The first one {Ii(z), Io(z)} corresponds to the expansion around 
the singularity z = 0 : 


h(z) = z -? (1 - z)-i 2 F, ttto. 1 + 5fS; j) , 

h(z) = z-i(i-z)- l i 2 F 1 ( j -y- t '+v-q t ' + m + ^;2). 

and they correspond to the following holomorphic fusion rules 


Hz) : 
h{z) : 


X $12 -> $ c 
X $12 ““y $ 


(3.65) 


(3.66) 


a 1 - 


The constraint a—a' ^ 2Z ensures that the two fields produced in the fusion <I> ai x < I ) i 2 
do not have conformal dimensions that differ by an integer. In turn this imposes 

2 In the Appendix, we comment on the case a — a' £ 2Z, for which logarithms appear in the 
conformal blocks. In the loop model at n = n this occurs for operators W mi oj with m£Z. 


18 











that the two conformal blocks I\(z) and I 2 {z) have different monodromies around 
z — 0: 

h(z) e-^hiz ), I 2 {z) -+ e~ iW I 2 (z ). (3.67) 

While I\(z) and I 2 {z) have Abelian monodromy around z — 0, they have a more 
complicated monodromy around z — 1. Another convenient basis { Ji(z), J- 2 (z)} 
corresponds to the expansion around the singularity z — 1, and has Abelian mon¬ 
odromy around it: 

Ji(z) -> e~ inb J\(z ), J 2 {z) -> e~ i7Tb ' J 2 (z) . (3.68) 

Their explicit expressions, which can be found in Appendix [Bj are readily obtained 
from I\{z) and I 2 (z) by exchanging z O (1 — z) and (a, a') (&,&')• These two 

bases are related through 

2 

h{z) = M w( a i b i d ) M z ) , (3.69) 

i=i 

where the explicit expression of the matrix M(a,b,c) is given in Appendix [Bj For 
the antiholomorphic side, we have similar expressions for the bases {A 1 (^), Z 2 (^)} 
and {Ji(z), J 2 (z)}, with cq —* dj. Since it obeys both PDEs, the full four-point 
function G(z, z) has to be of the form 

G (z, z) = Y^ x kk h(z) I k (z ). (3.70) 

k,k 


3.3.4 Consequences of locality 

Locality enforces the function G(z, z) to have a monodromy around z = 0 of the 
form 


G(z, z) -> e 2infl G(z, z). 


(3.71) 


Note that this phase factor is non-trivial when the fields < h ai ,a 1 and $12 are mutually 
semi-local. A typical exemple would be the spin and disorder operators for the Ising 
model, for which e 2l7rAl = —1. Since we consider the case when a — a' (f 2 Z and 
a — a! 2Z, the equality of phase factors e l7r ( a_a ) = e l7r ( a - a ') — e ^(a'-a) _ ^(a'-a!) 
can only be satisfied if X kk is a diagonal or anti-diagonal matrix. Since the change 
cxi —* a\ (resp. di —> d)) is simply a reparameterisation of the field $ aiiS 1 which 
exchanges the roles of h(z) and I 2 (z) (resp. l\ (z) and J 2 (z)), we can assume without 
any loss of generality a diagonal decomposition 


G(z,z) = X 1 I 1 (z)I 1 (z) + X 2 I 2 (z)I 2 (z ), 


(3.72) 


where we have renamed X u to X\ and A 22 to X 2 . These two constants are closely 
related to the OPE structure constants: 


Ad 

A 2 


A C ^$ 12 , < f> ai!Q , 1 , 4 > Ql _^h )(5i _^h^ C ^$ 0 , 3 ^ 3 , $ 0 ^ 4 , 04 , , 

A C ($ 12 , C ($a 3 ,a 3 , $ 04 , 04 , $ Ql + ^,a 1 + ^h) , 
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(3.73) 


where A is an unknown constant. 

It follows that the monodromy between $12 and an arbitrary held ^ h j l can only 
be e 2j7rAl = ±1. Indeed, the monodromy around z = 0 is 

G(z, z ) -G e i ^~ a) X l h{z)h{z) + e i ^'- a ' ) X 2 I 2 (z)I 2 (z ), (3.74) 


and generically we expect both Ad and X 2 to be non zero. This requires 2/i = 
a — a = a' — a' mod 2, which boils down tc0 


a — a G Z, 


(3.75) 


and therefore fi G Z/2. We prove in the Appendix that this condition has to be 
obeyed by all primary fields. Let us stress out the meaning of this relation. The 
decoupling of the null vectors of $12 puts strong constraints on the spectrum of the 
CFT : the conformal dimensions of the primary fields $ h ^ m ust be of the form 





Zrv 

h = a(a — 2ao), 

h = f3(/3 - 2a 0 ), 

with 

) CM 

1 

w 

1 

Q 


(3.76) 


While this holds trivially for spinless fields, for which we can choose a — /3, this is 
a severe constraint for non scalar fields. The full spectra of both the TL and O (n) 
loop models satisfy this condition, although $12 does not belong to the spectrum of 
O(n). This is perhaps not so surprising. Indeed, the O (n) loop model contains the 
field $ 13 , and the corresponding differential equation naively yields the constraint : 
h(a+a + ) — h(a+a + ) = h(a) — h(a) = h(a—a + ) — h(d—a + ) mod 1 , which boils down 
to a — a G Z. More generally, we expect the differential equation corresponding to 

to be consistent with locality as long as a —a is integer. From now on, we assume 
that all fields obey this constraint, namely (a — a) G Zck_/2 [or ( a' — a) G Za_/2]. 

As a side remark, if a CFT contains both $12 and $21 (and if their null vectors 
decouple), fields in the spectrum must belong to one of the two following families: 

• Scalar fields a = a, with a arbitrary. This is the case for all fields in Liouvillc. 

• Spinful fields a = a m , e , ot = a- m , e , with e ,m G Z/2, with spin s = —me G 
Z/4. 

While a continuum of scalar fields is possible, only a discrete set of spinful fields 
is allowed in that case. Let us stress that the spin fields of the loop models under 
consideration do not fall into this classification, since these models contain fields 
with m = —m 7 ^ 0 and e ^ Z/2. This is allowed because $21 is absent from the 
spectrum of the loop models. 

3 This comes from our choice of working with a diagonal matrix X. If one works with an anti¬ 
diagonal X instead the conditions becomes = a — a' = a 1 — a mod 2, or equivalently a' — a £ Z. 
But as we already mentioned this is a pure matter of convention. 
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3.3.5 Solution of the monodromy problem 


Let us go back to the four-point function G(z,z) of (3.72). Repeating the same 
arguments around z — 1 yields a decomposition on the (J^z), J 2 (^)} basis : 

G(z, z) = Y 1 Mz)Mz) + Y 2 Mz)J 2 (z) . (3.77) 

Compatibility with the change of bases yields two consistency relations : 


M\ i M\ 2 Xi + M 21 M 22 X 2 — 0 , 
M\ 2 M\ i X\ + M 22 M 2 \ X 2 = 0 . 


(3.78) 


For the the two linear equations ( |3.78[ ) to have a non-zero solution (Xi,X 2 ) (and 

(3.79) 


therefore a non-vanishing four-point function), we must have 

MuM 22 M\ i M 22 


det 


M\\M\ 2 M 2 \M 22 
Mi 2 Mi i M 22 M 2 i 


= 0 




M] 2 M 2 1 M\ 2 M 2 \ 


or more explicitly 


sin 7r ( p+a 2 fc sin 7r 


(p+b—a—d\ 

V 2 ) 


sin 7r 


(p+a—b—d\ 

V 2 J 


S1117T 


/p+b—a—d\ 

V 2 J 


SU17T 


/ p+d—a—b \ 

V 2 / 


sin 7T ( 3 P~ d ~ a ~ fc ) sin7r (P±W) siii7 r( 3 P -J-a-5 ) 


(3.80) 


where a —a, b — b and d — d are all integers. This constraint is manifestly symmetric 
under exchange of a and b, but not under exchange of a and d. This is because our 
analysis has singled out the relation between the conformal blocks around z = 0 and 
z — 1. Repeating the same steps around z — oo, we would end up with two more 
relations: one symmetric under a <(-)■ d and another symmetric under b <(-)■ d. When 
a, b and d are generic, all these constraints boil down to 


(a + b + d) = (a + b + d) mod 2 . 

The solution can then be written as 

M 2 \ M 22 M 22 M 2 \ 


Xi 

X 2 

or in a more symmetric fashion 


w 

AT 


M11M12 M12M11 


M 22 M 2 i M 22 M 21 

M11 M 12 M11M12 


(3.81) 


(3.82) 


(3.83) 


This leads to 


A, 

A 2 


7(^) 


7 


d—a— 


^)t( 


d—a—b+p 


) 


X (oij —> Uj) . 


(3.84) 


7 (^) 7 ( d ~ Q, 2 b ' +P ) 7 ( d ~ a ' 2 ~ 6+P ) 

This expression is invariant under b O 6', d O d' or b d, as it should be. 
In particular, taking b = b = —p and (d, d) = (a, a) we recover equation (3.57), 
obtained in the framework of the Coulomb gas. 
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3.3.6 OPE coefficients 


Using (3.57), (3.73) and (3.84), we can extract a ratio of OPE coefficients, and if we 


denote <f> 2 = $ Q3 ,q 3 and $3 = ®a 4 ,a 4 , we get 



(3.85) 


This formula immediately applies to the loop model CFT, since the latter satisfies 


the assumptions made in the beginning of this section. Like in (3.57), although it 
is valid in the presence of electromagnetic operators W me , it only gives a recursion 
relation on the electric charge. 


Equation (3.85) is a functional relation satisfied by the OPE coefficients. One can 
check that the RHS of ( |3.85| ) is what one would get if one replaced C(... < P Qii a i ...) 
by \/Cc{- ■ ■ V ai ... )Cc(- ■ ■ Vai ■ ■ ■ )• However, unlike the Liouvillc theory where two 
functional relations (coming from the degeneracy of <f> 12 and <f> 21 ) determine uniquely 
the solution Cc, in our situation the field $ 2 i is absent from the spectrum, and we 
only have one functional relation, which only determines the OPE coefficients up to 
a periodic factor. 


3.3.7 Application to the loop model 


When specialised to the TL loop model, (3.85) yields, in particular, the ratio 

1/4 


7(1 -pb(-i - p) 


C(W 10 ,W 10 ,W 12 ) _ 

Ciyv 10 , Wio, Wio) Lt(— 1 + 2p) 7(1 + 2 p ) 7 (—1 + p) 7 ( 1 + p) _ 


x 


7 3 (p/2) 7 (-l + p/2) 
_7 3 (-p/2) 7 (-l - p/2) 


1/2 


(3.86) 


which is checked numerically in Sec. 3 A Besides the precise expression of the ratio, 
it is important to notice the non-vanishing of the numerator: 


C(Wio,Wi 0 ,Wi 2 ) 7^0 


for 0 < n < 2 


(3.87) 


Indeed, since Wi 2 has conformal dimension h \2 in the holomorphic sector, one could 
have expected that it obeys fusion rules of the form Wi 2 x V a) a —>■ with 

(3 = a ± aq_/2, whereas this rule is clearly violated by C^Wio, Wio, Wi 2 ). As we 
shall see in Sec. [4j this is due to the presence of a logarithmic partner for Wio- 


3.4 Numerical checks 

Using a numerical algorithm based on transfer-matrix diagonalisation (see Ap¬ 
pendix 0. we have computed numerically the OPE coefficients in the TL loop 
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model at the point n = n, for 0 < n < 2. We can then compare these data to 


our predictions (3.55), (3.56), (3.57) and , (3.85), in particular in the presence of 


watermelon operators W me . 

3.4.1 The OPE coefficient C(Wio, Wio, V13) 

This OPE coefficient involves two spinless operators together with a degenerate 


operator of the form Vi^fc-i, and thus one expects it to be given by (3.25): 


C(W m , W w , V13) = J fihbh -Mlh ^ 7 ( p + l) 7 (p - 1) . ( 3 , 88 ) 

The agreement with numerical data on cylinders of circumference L — 10 sites to 
L — 18 sites is good, as shown in Fig. [5j 


3.4.2 The OPE coefficient C(Wn, Wn, V13) 


In this case we use the prediction of (3.55) with k = 2: 


c(v„ A , v„, a , v 13 ) = J2i(dlh_HdlP_JW 7 (i + o)7 (i + a')T(i + 0 ) 7(1 + o'), 


7(2p) 


(3.89) 


where 


a = 2a + a:, a = 2a + a:, a' = 2a + (2ao — a), a' = 2a + (2a 0 — a) • (3.90) 


The expression (3.89) is ill-defined when a —» an and a; —> and hence we need 

a regularisation procedure to extract C(Wn,Wn, V13). Indeed, if we set a = an + e 


and a = a_i,i + e with e, e —)■ 0, we get 7(1 + 0 ) 7(1 + a) ~ e/e. The condition (3.75) 
imposed by locality of the four-point functions yields a unique possible relation e = e, 
which finally leads to 


C{W U ,W n,V 13 ) = 


w 3 (p)p( 2 — 3p)y(l + 2p) 


7 ( 2 p) 


(3.91) 


The agreement with numerical data is good, as shown in Fig. [ 6 j The numerical data 
show an apparent divergence at n — 1 , but the height of the corresponding peak 
decays as L grows. We interpret this peak as a numerical artefact arising from a 
level crossing in the zero-leg sector with zero momentum, between the states |Vi 3 ) 
and T_ 2 -h— 210 ). 


3.4.3 The ratio C(V 12 , W n , Wi 2 )/C(Vi 2 , W n , W ; 


10J 


This ratio is predicted by (3.57), using the same regularisation procedure for Wn 
as above: 


C(V 12 ,W n,W 12 ) 
C'(V 12 ,Wi 1 ,W 10 ) 


7(-l + 2p)y(l + 2p)y(l - p)y(—1 - p) 


1/4 


7(“1 + P)7(1 + P) 


(3.92) 
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Figure 5: Numerical data for the OPE coefficient C'(Wio, Wio, V 13 ), compared to 
the expression (3.88). 


The agreement with numerical data is good, as shown in Fig. [TJ In this case, 
although h \2 + /i- 1,2 < 2/iio + 4 on the whole interval 0 < n < 2, in finite size some 
level crossing occursat n = n* between the lattice representatives of | W 12 ) and some 
fourth-level descendents of |Wio). This affects the numerical data close to n — 0, 
but the crossing location n* —> 0 as the system size grows. 


3.4.4 


The ratio C(W 10 , W 10 , W 12 )/C(W 10 , W 10 , W 10 ) 


This ratio is predicted by (3.86). The agreement with numerical 
shown in Fig. [ 8 j although the data still suffer from the finite-size 
| W 12 ) with the descendents of |Wi 0 ) discussed above. 


data is good, as 
level crossing of 


4 Indecomposability in the bulk theory of generic 
loop models 

4.1 The periodic Temperley-Lieb algebra 

Following the discussion in [13] , let us first describe the diagrams defining the alge¬ 
bra, and then give the algebraic rules between generators. We consider an annulus 
with L marked points on the inner boundary and L marked points on the outer 
boundary (L is an even integer). Actually, we shall depict the annulus as a rect¬ 
angle with periodic conditions across the vertical sides, and draw the inner (resp. 
outer) boundary at the bottom (resp. the top) of the rectangle: see Fig. [9j 
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Figure 6 : Numerical data for the OPE coefficient C'(Wn, Wn, V 13 ), compared to 
the expression (3.91). 



Figure 7: Numerical data for the ratio C(Vi 2 , Wn, Wi 2 )/C'(Vi 2 , Wn, W 10 ), 
pared to the expression (3.92). 


com- 


25 


























o 

g 

o 

g 

(M 

g 

o 

g 

o 

g 



Figure 8: Numerical data for the ratio C'(Wio, Wio, Wi 2 )/C'(Wio, Wio, Wio), com¬ 
pared to the expression (3.86). 



Figure 9: Two distinct diagrams of PTLl = 4 (see Dehnitions[fJ[2]and[3]). The top and 
bottom boundaries of the annulus are drawn as full lines, whereas periodic boundary 
conditions are assumed across the dotted lines. The parity of the left diagram is 
e — 1, and the parity of the right diagram is e = — 1 (see Definition [4]) . 
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Definition 1. A diagram w is a set of lines inside the annulus, connecting the 2 L 
marked points with no intersections. The arches (resp. the legs) of w are the lines 
connecting two points on the same boundary (resp. opposite boundaries). 

We can consider diagrams as basis elements of an associative algebra: 

Definition 2. The periodic Temperley-Lieb algebra PTL l is the algebra generated by 
the diagrams on an annulus with 2 L marked points, where any two diagrams (wi, W 2 ) 
are identified (up to a multiplicative factor) according to the following rules: 

• if u >2 can be obtained from w\ by a continuous deformation of the lines, then 
wi = w 2 ; 

• if w 2 can be obtained from W\ by a rotation of angle 2i r of one of the boundaries, 
then wi = w 2 ; 

• if w 2 can be obtained from W\ by removing a closed, contractible loop, then 
wi = n w 2 ; 

• ifw 2 can be obtained from w\ by removing a closed, non-contractible loop, then 
wi — n w 2 ; 

and the product w±w 2 of two diagrams is obtained by gluing the bottom (or inner) 
boundary of w\ onto the top (or outer) boundary of w 2 . 

Equivalently, we can formulate the above rules in terms of generators, which act 
by convention from bottom to top (or inner to outer) boundary: 

Definition 3. The periodic Temperley-Lieb algebra PTL^ is the algebra generated 
by ei,... eL and the right-shift operator r, subject to the relations for 1 < j, k < L: 

e j e j±i e j = e j ) 

e| = n ej , 

eje k = e k ej if 1 < \j - k\ < L - 1, 

r L = l, (4.1) 

T€jT ej- )-i , 

e\e 2 ... e k -\ = r 2 e ^i, 

(e 2 e 4 • • • e L )(e l e 3 ... e L -i) = n r(eie 3 ... e L _i), 


where e^+i = ei. 


Finally, let us define a notion of parity: 

Definition 4. The parity of a diagram w is e[w) = +1 (resp. e(w ) = —1 ) if the 
lines of w intersect the vertical boundary an even (resp. odd) total number of times. 

With this definition, we see readily that e^ii/^) = e{wi)e{w 2 ). 


27 



Wn 





1 2 3 4 

W 2 : 

12 3 4 





4 1 2 


1 2 3 


Figure 10: The basis states of W 0l W\ for L = 4 and W 2 for L — 8. 

4.2 Representations of PTL^ 

4.2.1 The spaces Wt and the TL scalar product 

Under the action of the above algebra, the bottom arches are left unchanged, while 
the number of legs never increases. This suggests to define the action of PTL^ on 
the “top-half” of the diagram only. For this, we now consider a circle with L marked 
points, with a outline connecting the origin to the exterior of the circle. Equivalently, 
we can place the points on a horizontal segment, with periodic conditions across a 
vertical cutline. 

Definition 5. A state a (also denoted \ a)) is a set of non-intersecting lines drawn 
in the lower-half plane, with some lines (arches) connecting (L — 2d) marked points 
among themselves, and 2d vertical labelled lines (legs) attached to the remaining 
marked points, so that the labels form a cyclic permutation of (1, 2,... 2d), and a leg 
attached to the point j e {1,..., L} carries a label with the same parity as j. 

In particular, the top-half of a diagram w, consisting in the arches on the top 
boundary of w together with the legs of w, defines a state. If we fix the number of 
legs, we get a representation of PTL l- 

Definition 6. The representation Wt is defined as the geometric action o/PTLl 
on the vector space generated by the states with 2d legs, where the contraction of any 
two legs gives zero. 

Definition 7. The space W = Wo © W\ © W 2 © ... is equipped with a bilinear, 
symmetric form (• , ■), defined as follows. Let a and b be two states with d a and db 
legs, respectively. Consider the set of lines g(a,b ) obtained by gluing the reflection 
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of a around a horizontal axis with b, and set 


{ 0 if any two legs of a are joined in g(a, b), 

0 if any two legs of b are joined in g(a, b), (4.2) 

n N c (a,b) ^N nc (a,b) otherwise, 

where N c (a,b) (resp. N nc (a,b)) is the number of contractible (resp. non-contractible) 
loops in g(a,b). In particular, if I a ^ £ b then (a,b) = 0 automatically. 

Proposition 1. The generators Cj are self-adjoint for the bilinear form (■ , ■) on 
the space W: 

V(a, b) G W 2 , (ejd, b) = (a, ejb). (4.3) 


4.2.2 The map ip 

Let us now describe a linear map whose properties will allow us to relate the spaces 
W 0 and IT,. 

Definition 8. Let ip be the linear map defined as: 



r it, -g w 0 




(p : < 

p-V-T.) - w +1 (i 

j...k 

... u ...)-| 

j...k \ 

... u ...)) 

(4.4) 


j ...n. j ...A/ 

where U (resp. U ) means that the points j and k are connected by an arch which 
intersects the cutline an even (resp. odd) number of times. 

We see immediately that 

j J + 1 ('_ 1')! +1 j j +1 

e M--\ I •••> = V < n "")!■■■ u ■■■>■ ( 4 - 5 ) 

Hence, taking all the signs into account, ip is a morphism of algebras, up to a parity 
sign: 

Proposition 2. If n = n, then 

\/w , Va G Wi , w ip\a ) = e(w) (pw\a). (4.6) 

Proposition 3. If n = h, then 

Va G Wo , VftGWi, (a, tp(b)) = 0. (4.7) 

In particular, ip{b) has zero norm: {(p(b), (p{b)) = 0. 

Note that a state with no legs can be specified by a planar pairing of the L 
marked points, with each arch decorated by the parity of its number of intersections 
with the cutline. In the case n = n, one may define a representation which ignores 
the parities: 
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Definition 9. A zero-leg reduced state |a) is a planar pairing of the L marked points. 
When n = n, the representation W 0 is defined as the geometric action of PTL l on 
the vector space generated by the zero-leg reduced states. 


With some simple counting, one gets the dimensions of the above representations: 

L\ 


dim Wo = 


L 

L/2-l ) ' 


dim Wn = 


(L/2)!(L/2 + l)! 

Proposition 4. If n = n, then p is injective, and 

W 0 — W 0 © Im (p) , Im (p) = W 1 . 


(4.8) 


(4.9) 


Proof. Consider a state a G W 0 - If a has no arch with odd parity, then a G W 0 . 
Otherwise, we consider the most external arch of a with odd parity, and denote by 
j < k the ends of this arch. We define 

\a) = \..M...), \a') = \... J U...), \b) = 2(-iy | (4-10) 


From (4.4), we have 

\a) = \a') + <p\b), (4.11) 

where a has one less arch with odd parity than a. By iterating this process, we can 
write a as 

l a ) = I °k) + p\bk) i d'k G W 0 , bk G W\. 


Hence, we see that 

Wo C Wo + Im (ip) 


dim W 0 < dim W 0 + rk(^) . 


But, from (4.8), we have 


dim W 0 = dim W 0 + dim W 1 , 

which leads to rk(<^) > dim IFi, and proves that p is injective. 


(4.12) 

(4.13) 

(4.14) 

□ 


4.2.3 Degeneracies of H between the Wo and W\ sectors 

Let us consider, for n = n, the Hamiltonian associated to the TL model: 

L 

H = (4.15) 

3 =1 

Proposition 5. To any eigenstate \ W) G Hj of H, one can associate a degenerate 
eigenstate |y) = G Wo, and |y) is orthogonal to W. In particular, (x,x) = 0- 

Proof. Consider an eigenstate |W) G Hj, with energy E\ 

H\W) — E\W). (4.16) 
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From Prop. [4j the state |x) = <£>|W) G Wo is non-zero. Moreover, since every term 
in H has even parity, we have from Prop. [2j 


H |x) = H<p |W) = pif |W> = S|x), (4-17) 


and hence |x) is also an eigenstate of H, with eigenvalue E. Note that, from Prop. [3j 
|x) is orthogonal to the space W 0 , and therefore it is orthogonal to every state in 
W: we say that |x) is a decoupling state. '0 


From the expression (2.18) of conformal dimensions, which determine the en¬ 
ergies in the continuum limit, we see that, for generic n, the above are the only 
possible degeneracies between two sectors Wa and Wy with t / f. Note that the 
situation is very different when g is rational (see e.g. iei)- 


4.2.4 Mixed representation 

Some of the correlation functions in the loop model cannot be expressed in terms of 
scalar products in W, but instead they impose the use of a representation in which 
the contraction of two legs is allowed. 

Definition 10. The representation W <i is defined as the geometric action o/PTL^ 
on the space Wo © W\, with the following coefficients for contractions: 


3 3 +1 


3 3 +1 


-3 I’ ’ ’ 
1 


...)=e j \... U ...), 


, 1...L, 


£l\ \ ■ ■ ■ \ ) — € l\ U ), 


(4.18) 


with ej = (—1) J+1 . 

Definition 11. The space W <i is equipped with a symmetric, bilinear form ((• , ■)), 
defined as 

((a, b )) = a(a, b ) a'(a, b ) n N ^ a ' b) n N ^ a ’ b) , (4.19) 

where the notations are the same as in the definition of (a, b), and 


cr(a, b) 


* 

(_X)Pai,+ta + l 

< 


if two legs of a, connected to the points j a < k a 
are joined in g(a,b ) by crossing the outline p a b times, 


1 otherwise, 


and cr'(a,b ) is defined similarly in terms of the legs ofb. 

Proposition 6. The generators ej are self-adjoint for the bilinear form ((• , •)) on 
the space W< \: 

V(a, b ) G W%, , «e,a , b)) = ((a , e,-6» . (4.20) 
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4.2.5 Jordan blocks 


Let us set n — n. For any pair of degenerate eigenstates |W) and |x) = (p\ W) 
discussed in Prop. [5j we shall argue that H has a Jordan block in the representation 
W< i, with matrix: 

( 0 E ) • < 4 ' 21) * 

Let us view |W) and |x) as vectors in the representation W<\. The action of H on 
|x) is unchanged, since |x) is in the zero-leg sector. For |W), we have to include the 
zero-leg terms arising after a contraction, which we denote |W C ): 

H\W) = E\W) + |W C ). (4.22) 


Suppose there exists a solution |W 7 ) of the linear system in Wq : (H — E)\W) = 
|x) — |W C ). We have checked numerically the existence of a solution^ in many cases, 
including all the pairs of degenerate states mentioned in this paper. Note that for 
any scalar A, the change |W 7 ) —> |W 7 ) + A|x) does not affect the above linear system, 
and thus the solution is not unique. We introduce the states 

lx» = lx) . |w» = |W) + I W) , (4.23) 


and we get 


H\x)) = E\ X )), H\W)) = E\W)) + \ X )). 


(4.24) 


In summary, given an eigenstate | W) with energy E in the two-leg representation 
W \, we have constructed a pair of states (|x)), |W))) in the representation W <i so 


that E[ has the Jordan form (4.21) in span(|x)), |W))), and |x)) is a state in the 
zero-leg sector of W<\ with vanishing norm: ((x, x)) = 0. 


4.3 Scaling limit 

Throughout this section, we set n = n. 


4.3.1 Identification of states 


The states {|Via,)} and {|W me )} defined by vertex operators in Sec. 2.3 can be 
identified as the scaling analogs of states in the representations {Wi}. The general 
idea is that the continuous analog of PTL^ is made of two copies of the Virasoro 
algebra: Vir © Vir, and that the discrete quantum numbers of E[ survive in the 
scaling limit. 

For clarity, we denote |$)l and (T)) i a state in Wi and W <i respectively, and 
1$)) the corresponding state in the CFT Hilbert space. We use a similar notation 
for scalar products. 


4 Note that the signs e 7 in the action of H on W<i are crucial, since the momenta of |W) 

and |x) = y|W) differ by n. If instead we had set ej = 1 (and a(a, b ) = a '{a, b) = 1 in the 

corresponding scalar product so that H is still Hermitian), then H would have a diagonal form in 

the block (|x», |W))). 
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• The state |Via;)l, with k > 1, belongs to Wo C Wo, and it has finite norm 
(which is set to one by convention), since the scalar product (•, •) is non¬ 
degenerate in Wo for generic n. 

• The state |Vi^)x,, with k < 0, belongs to Im<^ C W 0 , and it totally decouples: 

Vk< 0, Va G W 0 , (V lfc ,a> L = 0. (4.25) 

• The state |W me )L, with m > 1 and e G Z, belongs to W m . Recall that the 
allowed electric charges are of the form e = p/m, with p G Z, which may be 
interpreted as follows. Consider the cyclic translation R of the leg labels by 
two units: 


R\... | ... | ... | ...) = |... | ... | ... | ...}, ( 4 . 26 ) 

x\ X2 X 3 Xl+2 X2+2 x 3 +2 

where the labels are considered modulo 2m. The action of H on W m commutes 
with R , and we can label the eigenstates of H according to the eigenvalues 
r q = exp(2f7rg/m) of R, with q — 0,..., m — 1. Then the state |W me )L with 
e = p/m is in the sector r q such that p = q mod m. Again, since the scalar 
product (•, •) is non-degenerate in W m , one can set (W me , W me )L = 1- 

• For m = 1 and k > 1, Wi k is degenerate with the null descendant of Vifc 
(under the action of Vir) at level k, which we denote Xik■ In the representation 
Wo © W\, the two corresponding states (|xife)n, |Wi k)L) are related by 

\Xik)L = v\Wik)L, ( 4 . 27 ) 


and they are independent, orthogonal eigenstates. In contrast, in W< i, the 
corresponding states form a Jordan block, and we use a disctinct notation 
(|Xifc))i) |! / Vifc))n) I° r them. The Jordan block (4.21) becomes, in the scaling 
limit: 

hik + hi : -k 1 

0 hik + h\-k 


Lq-^-Lo — 


(4.28) 


There is a similar relation between Wi and xik, the null descendant of Vik 
under the action of Vir. 


Moreover, as we have argued above, any eigenstate in W\ produces a 2 x 2 
Jordan block in W<\. In particular, since <p commutes with the action of PTL^ 
(which generates Vir © Vir in the scaling limit), we expect that in the scaling 
limit 


<pL_ ni .. ■ L _ np | Wife) — T-m • • ■ L- nv ip\Wik) — L_ ni ... L_ Up \xik) (4.29) 

and thus L_ ni ... L_ np \Wik )) and L_ ni ... L_ n \xikj) w ih also form a Jordan 
block. 

• Finally, let us comment on the case of Wio- Here, the vertex operator Vio is a 
null state, and in W <i we have a Jordan block for the pair (|Vio))l, |Wio))l). 
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4.3.2 Hilbert space 

A crucial question is then to decide which representations should be included in 
the Hilbert space, so that the scaling limit results in a consistent CFT, i.e. with a 
consistent operator algebra. Indeed, several unequivalent representations correspond 
to the same spectrum in the continuum limit. 

Let us first examine the choice W = Wq © W\ © W 2 ... as the Hilbert space, with 
the associated scalar product (•,•). I 11 this representation, the zero-norm vectors 
span the subspace Im ( 99 ) C W 0 , and they actually decouple from the whole vector 
space W. Once they have been removed, we are left with W = Wo © W\ © W2 • • •, 
where the scalar product (•, •) is non-degenerate for generic n. Hence, in the cor¬ 
responding CFT in the scaling limit, the 2- and 3-point functions will have the 
standard form, and the ordinary differential equations arising from null-vector con¬ 
ditions will hold. These two facts are sufficient to ensure fusion rules of the form 

x 3*1,71 ~t $ li | n _ m | +1 + • • • + , (4.30) 

separately in the holomorphic and anti-holomorphic sectors. However, the results 
exposed in the end of Sec. |3.3| show that these rules are not obeyed, e.g. the fusion 
W 12 x W 10 —>• W 10 is allowed. 

This contradiction is resolved if we use instead: 


n = w< x © w 2 © w 3 ©... 


(4.31) 


with the associated scalar product ((•,•)). In simple terms, it means we take the 
same basis of states as in W, but we allow the mixing of some states from the zero- 
and two-leg sectors, both in the action of the algebra and in the scalar product. 

In the following, we shall assume that the representation (4.31) is the one that 
produces a consistent operator algebra in the scaling limit. 


4.3.3 Decoupling states 

We are now ready to discuss which of the zero-norm states actually decouple from the 
whole Hilbert space % under the scalar product ((•, •)). Let us denote by A ik e Vir 
the combination of generators which produces a null descendant at level k in the 
module of V\ kl and A lk the analogous object in Vir: 

Ixifc) = A lk \V lk ), \xik) = A lk \V lk ). (4.32) 

As explained above, the zero-norm state |xifc) (resp. |xifc)) forms a Jordan block 
with |Wi_fc) (resp. (Wife)). However, this fact does not tell us anything about the 
scalar products (xik, Wi - k ) = {xik, Wife), and at present we are not in a position to 
make any statement on the decoupling of the vectors \xik) and |xi*) ■ We did check 
numerically that the scalar products ((xifc, W’lkfjL = ((Xifc, Wi ,-kf)L for k — 0,1,2 
(setting xio — X 10 = V10) do not vanish at finite size L — 2,..., 12, but considerable 
additional work is needed to conclude about the scaling limit of these scalar products. 
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The operator Wu- (resp. Wi _&) also admits a null descendant AikWik (resp. 
AifcWi.-fc), since it has conformal dimensions (h\k, hi-k) [resp. (hi _fc, hi,fc)]. Acting 
with ip in the scaling limit, we get 

Wife)) = </5(A lfc |)Vi _fc)) = A lfe A lfe |Vife), (4.33) 

and thus A lfc |Wifc) = A lfc |Wi,_fc), since ip is injective. Indeed, in the lattice model, 
we find a unique state in W\ with energy ~ ^ x 2h 1 _/ c and zero momentum, de¬ 
generate with a unique state in W 0 . In the scaling limit, the scalar product between 
the corresponding states is ((AikAikVik, Au-Wifc)) = 0, because A\ k A lk \ Wife)) = 0 
by construction. Hence, the states Au,|Wifc)) = Aifc|Wi _&) and AikAik\V\kf) are 
totally decoupled, and we may set them to zero in the Hilbert space of the CFT: 

AifejlVifc)) = A lk \ Wi _fc)) = 0 , and A lk Ai k \V\kf) = 0. (4.34) 

Note that, in the case k — 1, the lattice analogs of Hi,±i are two of the “discrete 
parafermions” [2T, 22] of the TL or 0( n) model, i.e. lattice operators satisfying a 
discrete version of the Cauchy-Riemann equations <9Hfi,i = 0 and dW\-\ = 0. 


4.4 OPEs in the presence of Jordan blocks 

4.4.1 Limiting procedure 


The Jordan block structures discussed in previous sections appear at the point n = n, 
for generic n. We can thus follow a treatment of the singular OPEs similar to the 
c —> 0 limit of [71 [23], except that in our case we have a two-parameter family of 
CFTs (labelled by g and e 0 in the notations of Sec. [ 2 j) : we shall fix the coupling 
constant g , and let the background charge eo tend to (1 — g). 

More specifically, we have seen that for general parameters (n, n), we get for the 


TL loop model the central charge and dimensions (2.18) of the form: 


fie 2 I / p \ 2 p 2 

c(g,e 0 ) = l--k, hme (g,e 0 ) = l [rn^-— J (4.35) 

where m — 0, e e e 0 + Z for electric operators Vo e , and m e {1, 2, 3 ... }, e G Z/m 
for electro-magnetic operators W me , and 


n = —2cos7rg, n — 2 cos 7 re 0 . (4.36) 

In contrast, the dimensions of the Kac table for central charge c(g , eo) are of the 
form 

hlT ] ( 9, e o) = (ra_ + sa+f - j -, (4.37) 

4 g 

with (r, s) integers, and 


a + + cm =-, ct + a_ = —1. 

V9 


(4.38) 


The general idea which we develop in the following is to consider an OPE of two 
operators, such that one of the terms in the expansion is a null state belonging to 
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a Jordan block of L$ + Lq for n = n, and analyse the behaviour of OPE coefficients 
in the limit n —> n. 

In particular, the following dimensions of vertex operators coincide with the first 
row of the Kac table in this limit: 

ho,e 0 +k -t h[ K k+i t hi }k+ i —>■ h [ K k } +1 . (4.39) 


4.4.2 OPE coefficients 

Consider the OPE between two primary fields $1 and $ 2 , with dimensions (hi, hi) 
and (h 2 ,h 2 ): 

$1 (z, z)$ 2 (0) = $ 2 , $p) z - hl - h2+hp z-~ hl -~ h2+ ~ hp [<h p (0) + ...], (4.40) 

p 

where the sum is over all the primary fields appearing in the fusion $1 x <h 2 , and 
the ... denote the contributions from the descendants of <L p . In the following, we 
shall need the precise expression of the OPE coefficients of some of these descendants. 
By the standard methods of CFT, we obtain 

C(* 1 , $ 2 , L _ kl ... L_ km %) = 1 , $ 2 , *,) x , (4.41) 

where the /Jf^'satisfy a set of linear recursion identities. At level one we have 

2h p fii2 ^ = hi — h 2 + hp . (4.42) 

At level two, the coefficients of L 2 _ i<f> p and L_ 2 < h p are determined by the system 

2 ( 2 h p + 1 )/Jf 2 ' + 3/3^2 ■* = (hi — h 2 + h p + 1 )/Jf 2 ^ 

6h p /3f 2 ’ + (4h p + c/2)0^ ^ = ( 2 hi — h 2 + h p ), 

whose solution may be written as 

op( 2 ) -P(hi, h 2 , h p ) 

1 12 ~ 16h 2 p + 2(c - 5)h p + c ’ 

aP ( 1 , 1 ) _ 2 ( 2 hi — h 2 + h p ) — ( 8 h p + c)/?^ ^ 

P 12 — Tbh > 


(4.43) 


(4.44) 


where 


P(hi, h 2 , h p ) — —3(hi — h 2 + h p )(hi — h 2 + h p + l) + 2(2h p + l)(2hi— h 2 + h p ). (4.45) 

4.4.3 Null states in the OPE 

As we have seen in previous sections, at n = n, the null descendant Xik at level k 
of the primary operator V\ k may survive in the theory, forming a Jordan block in 
L 0 + L 0 with its logarithmic partner Wi Let us analyse in detail the case of Vi 2 , 
which is the energy operator in the TL loop model. 
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Consider two primary fields <f>i and $ 2 , such that their fusion contains the terms 


$1 x $2 — y Vo,e 0 +l + VV1-2 + ■ • • 


(4.46) 


where the ... denote primary fields with distinct dimensions. In the limit n —> n, 
we have 


h — h 0:£o+ i —> h[ 2 , h! — h\2 -> h[ 2 ^ ■ 

The contributions from level-two descendants of 14 = Vo, eo +i in 
be organised as 


(4.47) 

the OPE flOo] ) 


can 


P(h 1 ,h 2 ,h ) 


16h 2 + 2(c — 5)h + c 


L_n 


8 h + c 
12 h 


L -1 ) 14 + 


2 h\ — h 2 T h 
6 h 


L\V h 


(4.48) 


The numerator P(hi,h 2 ,h) of the first term is a polynomial which vanishes if h = 
h\ K p and (hi, h 2 ) obey the standard fusion rule with $12 (expressed in vertex charges 
hj = a 2 — 2a 0 aj) a 2 = ai±a+/2, whereas the denominator has roots h\P and h^P. 

If the dimensions hi and h 2 obey the standard fusion rule with $12 in the limit 
n —>■ n, then all the coefficients in (4.48) remain finite. As n — )■ n, the first term 
converges to the null state 


{L - 2 — gL 2 _ 1 ) V 12 = X 12 • 


(4.49) 


On the other hand, more interestingly, we shall now show that the presence of 
the operator Wi _ 2 allows the fusion x $ 2 V 12 , even if a 2 4 a 1 ±a+/2. In this 

case, we define the finite quantity 


Q 12 — 


-P(hi, h 2 , h) 
16 (h — h. 


(K)i 

21 t 


(4.50) 


and (4.48) takes the form 

Q 12 


h — h 


(K) 

12 


(L- 2 -gLl t )V„ + # LltVn, 


(4.51) 


where 4 denotes a finite coefficient. Hence, in the limit n —> n, the null state X 12 
appears in the OPE of dq x $2 with a diverging coefficient. This diverging term in 
the OPE (4.40) should be cancelled by a term coming from W4_2, which tells us 
that the latter operator should decompose as 


C(4 1 ,4> 2 ,V k ) , Q n 

VVl -2 — ~-T x 


C($l,$2,Wl,_ 2 ) ' h~h\ K) 


\— (L_ 2 — gL 2 _ 1)14 + k ^ 1 ,- 2 ] , (4-52) 


''12 


where wi- 2 has a finite limit as n —> n, and k is a vanishing coefficient. Note that 
wi _ 2 is not defined uniquely, since it we can add a term oc (L _ 2 — gL^^Vh to it. 
We expand the power factors which multiply the term (4.51) in the OPE: 


~-hi-h2+h+2--hi-h2+h _ 


z _ ftl _ fc2+ v+ 2 ~fc 1 -ft 2 +fc' x []_ + (/, _ h >} \ og (zz) + o (n - n)] . 


(4.53) 


37 
















Then, if we set the normalising factor in (4.52) to k — h — h', the OPE (4.40) takes 
the well-defined form 


$1 (z, z)® 2 { 0) = z - hl - h2+hl2+ ' 2 z- hl ~ tl2+hl2 C($ 1 , $ 2 , Wi _ 2 


_ ,— hi — h,2+hi2+2 ^—hi—h2+hi2^ 

x {wi - 2 ( 0 ) + [(L _2 - ^ 1 )Vo, eo +i](0) log(zz)} + 


(4.54) 


where the ... denote contributions from other vertex operators and other descen¬ 
dants, and the OPE coefficient is given by 


C(G> 1 , 4*2, 10 ,- 2 ) — C(<1?1) <& 2 ) V/J x Q 12 x 


h-h' 


h-h 


(K) 

12 


(4.55) 


Finally, combining the definition (4.52) of 1111,-2 and the well-known two-point 
functions for n^n: 


(V h {z, z)V h ( 0)) = \z\~ 4h , (Wi _ 2 (z, ^)W 1 ,_ 2 (0)> = |z|- 4 ' 1 ' 

(v fe ^,^)yv 1 ,_ 2 (o)) = 0, 

we obtain the two-point functions in the limit h —> n: 


(4.56) 


(*12(3, -)xi2(0)) = 0 , {wx- 2 {z, *)wi,- 2 ( 0 )) = 


9 - 2/3 log (zz) 

g2hi,—2 g2hi2 5 


(xMz, z)w l - 2 (0)) = 


p 


2h\ —2 ‘y2h\2 5 


Z 


(4.57) 


where 6 is arbitrary, and the indecomposability parameter f3 is given by 

(V,|(L 2 -^)(L_ 2 .^L 2 _i)|V ft ) 


(3 = lim 

eo-tl-g 


h - fl' 


where h = /io,e 0 +i and h' — hi 2 . A little algebra gives the final result 

0 = 4(1 -g). 


(4.58) 


(4.59) 


5 Conclusion 

In this work, we have shown how to adapt the conformal bootstrap approach to the 
case of operators with conformal spin, and applied it to compute some OPE coeffi¬ 
cients involving some mixed electric-magnetic operators in the O (n) loop model. It 
turns out that, in all the cases where this approach applies, the result is particularly 
simple: the OPE coefficient is the geometric mean of the timelike-Liouvillc DOZZ 
formula mm on the holomorphic and anti-holomorphic sectors. Our results are 
supported by economical, yet precise numerical calculations based on the transfer 
matrix of the loop model. In doing so, we used a numerical method which can be 
exploited more extensively to test many predictions on OPE coefficients. 
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Many OPE coefficients remain out of reach for this conformal bootstrap ap¬ 
proach, including some very natural observables from the geometric point of view, 
e.g. the probability that three points lie on the same loop. We have traced these diffi¬ 
culties to the non-diagonalisability of the dilatation operator (L 0 + L 0 ), and proposed 
a simplistic analysis of the representation theory of the periodic Temperley-Lieb al¬ 
gebra, along the lines of (9j HDJ HU H2| 1T3] . to account for this effect at the lattice 
level. Similarly to the representation theory of quantum groups, this problem is 
much simpler for generic values of the loop fugacity n = q + q _1 when q is not a root 
of unity, which is the case under study, but many questions still remain unsolved, 
the most important one being: do the null vectors descending from electric oper¬ 
ators |Vifc) become decoupled or not in the scaling limit? Answering this type of 
questions would require a substantially deeper analysis of the loop model for generic 
loop fugacity, which we leave for future work. 
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Appendix 

A Details of the Coulomb-gas approach 


The normalisation factor appearing in (3.7) is: 

k -2 


74(a:i, a 2, o?3, cut) — PJ 


3 =0 


r[(j + 1 )p\ r(i + a + j»r(i + c + jp) 

T(p) T[2 + a + c + (/c — 2 + j)p\ 


x 


p—k—1 

n 

j =0 


r[(j + 1 )p\ T(i + b + jp)T(i + d + j P ) 

r (p) r[2 + b + d+ (p - k - 1 + j)p\ 


(A.l) 


The matrix elements on the last row and last column of A [see (3.10)] are: 


k -2 


s(l + b + jp) 


p—k—1 

n 


s(l + d + jp) 


Apk s[2 + b + c+ (k - 2 + j)p\ ^ s[2 + a + d + (p - k - 1 + j)p\ 


, (A.2) 


p-2 k—2 

A kp =JJs[(j + l)p] 


s[l + b + (p - k + j)p\ 


j=0 / = o S [(J + !)P] s[6 + c + (p - 2 + j)p] 


x 


s[l + c+ (k - 1 +j)p] 


p—k—1 

n s[(i + l)p] s[b + c + (p - 3 + k + j)p\ ’ 


(A.3) 
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Figure 11: Integration contour defining: (a) the conformal block (b) the 

conformal block Fn{z). 


where s(x) = sin(7nr). The general 
useful to notice that 


' k—2 

A kp oc <{ + l)p] 

0=0 


p—k—1 

;(1 + c + jp) JJ s[(j + l)p]s(l + b + jp) 

3 =0 


(A.4) 


We shall give here the explicit calculation leading from the expansion coeffi¬ 
cients (3.23) to the expression (3.25) of the OPE constants. First, let us notice that 
the case k = 1 corresponds to C 2 (V a , V a , Vn) = 1, since Vn — 1- If we set a = ai p , 
and denote a p = 2a + ai p and a' p = 2a+(2ao — «ip), we get 


C 2 (v v v s s i\ a P, a 'p) 

G V '1,2k—1; Vlp, Vl p ) - -TpT -— 

s 1 (a p ,a' p ) 


_ TT 7[Q + l)p] TT l[(j + l)p] 7[-l + (p + 1 + j)p] 

f = 0 7[-l + 0 + 2 )p] f = o 7[O' + k )p\ 7 [ip - 1 - j)p] 


P yr 1 7[O' + l)p] 7[~1 + (p + k ~ 1 ~ j)p\ 
} = o 7 2 [(p — 1 — j)p] 


(A.5) 


Going back to general cc, we have 

o(p) ( A 

C(V.,V a ,V 1 ,w)C(V 1) M,V 1 „V„) = 1 ' , 

Sf ; (a, a') 

which gives, combined with (A.5), the result (3.25). 


(A.6) 
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B Hypergeometric conformal blocks 

In this section, we recall the different bases of solutions of the PDE 


gd\ 


1 

-1- 

Z Z 


d z + ( — 


h. 


+ /12 — hi 


z(z - 1) 


/(*-)= 0 . (B.l) 


Like in Sec. 3.3[ we use a Coulomb gas parametrisation of the conformal dimensions 
hi = a l (a i — 2a 0 ), and we introduce 


a — 2a + a \, b — 2 tt_|_ct 3 , cl — 2 tt_|_ct 4 , 
a' = 2a + a\ , b' = 2a + a' 3 , d! = 2a + a' A . 


(B.2) 


B.l Generic case : a — a! and b — b' 2Z 

Several bases of conformal blocks can be considered. In Sec. |3.3 we refer to two 
such bases. The first one {I\(z), ^(z)} has simple monodromy as the variable 0 
loops around 0 and corresponds to the expansion around the singularity z — 0 : 

h(z) = Z~S( 1 - z)-i 2 U d '-“- b +P -1 + z) . (B.3) 

h{z ) = z-‘i( 1 - z)-'i 2 F 1 ( d '-' 2 - t ' +p , d - a '? +p ; 1 + ^f-,z) , (B.4) 

while the second basis { J\(z), ^(z)} has simple monodromy around z — 1: 

Ji{z) = z~H 1 - z)~* 2 F 1 (^ ±£ , d ^±p- 1 + 1 - , (B.5) 

J 2 (z) = Z"T(l-z)-£ 2 F, ( d '- a ' 2 - b/+P , d - a '- 6, +^ ; ! + 1 - z) . (B.6) 


It is useful to know how these conformal blocks transform under the various repa- 
rameterisations O a': 

• Under the exchange a -H- a', we have h(z) <(->■ h(z), while J\{z) and J2(z) are 
left invariant^] 

• Likewise, under b o b', we have J\{z) -H- Jo(z), I\ and h(z ) being left invari¬ 
ant. 


• Finally, under d <(->■ d', the conformal blocks {/i(^), h(z)} and {Ji(z), Jo(z)} 
are invariant. 

The change of basis is given by the matrix M(a, b, d): 

2 

h(z) = ^2 M ke(a, b, d) J e (z) , (B.7) 

1= 1 


5 To see this, the relation 2^1(0, b\ c; z) = (1 — z) c a b 2 -Fife — a, c — b; c; z ) comes in handy. 
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where 


M (a, b, d) = 


r(i+^)r(^) 

p ^ d—a — b f -\-p ^ p ^ d f — a — 

i+^)r(i^) 


r(i+^) r (^) \ 

) r ^ i 


^ p ^ d—gt — b'+p ^ p ^ d' — a' — b / -\-p 


7 ^( d—a-b-\-p \ T-i / d / —a —6+p 

V 2 ) L \ 2 

r(i+syi)r(if£) 


(B.8) 


) r ^-a'-i,+,. ^ r ^ d'-a'-b+p ^ y 


As long as all four conformal blocks are well defined (i.e. a' — a ^ 2Z and b' — b ^ 2Z), 
this matrix is invertible, with inverse M _1 (a, 6, d) = M(b,a,d). 


B.2 Logarithmic case 

Suppose now the indicial equation around the singularity 2 = 0 is degenerate, namely 
a — a' € 2Z. Without any loss of generality we can assume a' — a = 2m where m is 
a positive integer. Since we have a' + a = 2p — 2, this means 


a = p — 1 — m, a' = p — 1 + m . 


In that case, the conformal block given by (B.3) is still a solution, 

h(z) = 2 1± ^( 1 - z)~* 2 F 1 ( d ~ b ^ n+1 , df-b+m+l . 1 + m . _ 


(B.9) 


(B.10) 


On the other hand, (p.4) is no longer well-defined, since its third argument is a 


non-positive integer. We can approach the degenerate case as 
a = p — 1 — m — e, a 1 = p — 1 + m + e, 


(B.ll) 


with small e. Let us denote {Ii <e (z), l 2 , e {z)} the corresponding solutions. The first 
conformal block encounters no issues as e —> 0 : 


lim J le (z) = h(z ). 

e—>-0 


(B.12) 


On the other hand, the second conformal block, as given by (B.4), has a first order 
pole 


h,e{z) ~ 


T(1 — m — e) fd—b — m-\-1 \ fd'— b — m + 1 


ml 


h(z ), (B.13) 


where (x) n = x{x + 1) ... (x + n — 1) is the Pochhammer symbol. For instance the 
second conformal block can be chosen as the limit 


h(z) = 


-l) m rn\ , 

-lim 

2 e-s >0 


p | d— a'— 6+p^ p ^d' — a'—b+p^ 


r(i + ¥) 




p | d—a'—b'+p J p / d!—a'—b' +p 


r(i + ¥) 


-M 1 ) 


(B.14) 


with a and a' given by (B.ll). Note that this expression still makes sense for 


b' — b E 2Z, since the r.h.s. has a well-defined limit as b' — b —> 2n , n G Z. Moreover 
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we have chosen the second conformal block to be invariant under b <(-)■ b'. This gives 
a solution of the form 


h(z)=h(z) log z + H(z) z {1 - m -' ,)/2 , 


(B-15) 


where H(z) is a regular function around z — 0, whose explicit expression we will 
not need. The monodromy around z = 0 is no longer diagonal: 


h{z)^e~ ma h(z ), 


I 2 (z) —> e 


2mli(z) + I 2 (z) 


(B.16) 


If b' — b ^ 2Z, the basis {Ji(z), J 2 {z)} is the same as for the generic case (B.5)-(B.6), 
and the change of basis becomes 


M(a, b, d ) = m\ 


( 


\ 


\ 


d— b^ + l+^^p^ d f — b f -\-\-\-m ^ 

( ^^mp^ cZ—b+1 —m ^p^^ —b+1 — 

2r(i+^) 


p/ d-b-\-l-\-m ^p^ d f ^ 


(-i) m r( d ~ b '+ 1 ~ m )r( 


d'-b' + l-m 


2r(i+^#) 


M 1 (a, b,d) = —- 


1 

m\ 


/ p ^ d—b f -\-l-\-m ^ p ^ ^ ~b f 

2r(^) 

p^ d—b-\-l-\-m ^ p ^ d/_ 


-p+l+m 
2 


— b+l+m 
2 


V 


2 r(^) 


(-l) r 

(-1)" 


i+V) \ 


p ^ d— b-\-l — m ^ p ^ d! — b-f-1—m ^ 


' d-b' + l-r. 


^) r ( d '- 6 ' + 1 - m ) 


C Constraints on the spectrum from locality 


In this appendix we explain how locality imposes the constraint : 

a — aGZ or a' — aGZ. (C.l) 

We first prove this in the generic case (a' — a ^ 2Z) before doing the log case 
(a' - a E 2Z). 


C.l Generic case 

Let us consider the four-point function 

G(z,z) = (t ft „ il (0)4 1 , 2 (j)4 1 , 2 (l)4 klJ - 11 (oo)). (C.2) 

Note that this four-point function cannot vanish due to the trivial fusion channel to 
the vacuum. To describe this four-point function we can choose b = b = 3p — 2 and 

d = a. If there are two fusion channels, then we must have a — a G Z or a' — a G Z as 

we explained in the main text. Suppose now there is a single fusion channel around 
z = 0. Using the freedom to change a —)• a', we can assume 

G(z,z)=X 1 I 1 (z)I 1 (z). (C.3) 
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Since there is a single fusion channel as z —> 0, there is also a single fusion channel as 
z —» l[^] There we know that this fusion has to be $ 1,2 X $12 —> 1, which corresponds 
to 


G(z,z) = . 


(C.4) 


Consistency of these two relations enforces M 12 = M 12 = 0, where M is given by 
(B.8) with b = 3p — 2 and d = a. This requires respectively a G Z and a € Z, and 
in particular a — a G Z is obeyed. 


C.2 Logarithmic case 

Let us now consider the logarithmic case, namely a' — a = 2m with m 6 N. This 
means a = p — 1 — m and a' = p — 1 + m. If we also have a' — d = 2m with rh 6 N, 
then a — d = rh — m € Z. Out of interest we mention the form of the four-point 
function compatible with locality: 


G(z,z ) 


X 1 / 1 (^)7i(z)+X 2 




(C.5) 


with monodromy factor e. m ( a ~ a ) as in the generic case. 

This only leaves the hybrid case to consider : a logarithmic holomorphic side, say 
a' — a = 2 m with m e N, coupled to a generic antiholomorphic behavior, a' — d (f 2Z. 
In this case locality demands a single fusion channel, which we can choose to be 

G(z,z) = X 1 I l (z)I 1 (z). (C.6) 


The same arguments as in the generic case apply, and we find that a and a have to 
be integers. This is impossible since a = p — 1 — m. Therefore the hybrid case is 
ruled out by locality. 


D Numerical determination of OPE coefficients 
from the transfer matrix 

In this section, we describe a simple algorithm to compute OPE coefficients numeri¬ 
cally using the eigenstates of the transfer matrix or Hamiltonian. Only the principle 
of the method is given here, and we will expose all the implementation details for 
its application to loop models in a future publication. 

Let (0i, 02, 03) be three quasi-primary operators of a OFT, and suppose we want 
to determine numerically the OPE coefficient 0(01,02,03). This coefficient can be 
expressed as the scalar product 

C(01,02,0 3 ) = <01102 |03) • (D- 1 ) 

Let us consider the transfer matrix O for a periodic system of circumference L , 
denote by |0j)z, the eigenstate of O corresponding to the state 0j) in the OFT, and 

6 From the matrix X we get the expression of the matrix Y = XPXM. These two matrices 
share the same rank, and the number of fusion channels around z = 0 and z =1 has to be the 
same. 
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suppose that \4>i)l and \<Pz)l are obtained numerically by a partial diagonalisation 
of tL. Moreover, suppose that one is able to construct a lattice operator (p[ L \r) 
which behaves as (p^ ~ N 2 L~ X2 02 hi the scaling limit, where N 2 is some non- 
universal normalisation constant. More precisely, there will exist a series of operators 
02, 02 ■ ■ ■ °f increasing scaling dimensions X 2 < x 2 < ... and some normalisation 
factors N 2 , N 2 ,... so that 

0^ } = N 2 L~ X2 02 + N' 2 L~ x ' 2 0' + ... (D.2) 


in the scaling limit. Clearly, with the standard normalisation of CFT states (0^ |0y) = 
1, and denoting |0) the ground state, one has 0 2 L ' ) |O )l ~ N 2 L~ X2 10 2 ). Then one 
obtains: 


n(0l|02 L) |03)L 

l(02|02 L) |O)l 


C(0!, 02, 0 3 ) + 0(1/L X 2- X2 ) , 


(D.3) 


where x 2 is the scaling dimension of the most relevant subleading operator in the 
expansion (D.2) contributing to l(0i|02 L ' > |03)l or l(02 |02 Z |O)l- 
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